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ABSTRACT 

'o\ : 

We have developed a computer model to calculate gamma ray burst (GRB) light 
t-h I curves and efficiencies from the interaction of a single, thin blast wave with clouds 

in the external medium. Large amplitude, short timescale variability occurs when 
1—5 ' the clouds have radii r <C R/T, where R is the mean distance of a cloud from the 

^f) • GRB source and T is the blast-wave Lorentz factor. Efficiencies ^ 10% require a 

large number of small clouds, each with sufficiently large column densities to extract 
fS| . most of the available blast-wave energy in the region of interaction. The number and 

duration of pulses in the simulated GRB light curves are compared with the respective 
I ! properties found in GRB light curves. If GRB sources are surrounded by clouds 

with such properties, then short timescale variability of GRBs can be obtained in the 
' external shock model. 

OO : 

On ■ 

Or 

6 ■ 
a: 

23 ' 1- Introduction 

• i-H ! The cosmological origin of GRBs has been established as a result of optical follow-up 

observations of fading X-ray counterparts to GRBs discovered with the Beppo-Sax mission (e.g., 
Costa et al. 1997; van Paradijs et al. 1997; Djorgovski et al. 1997). The total energy released in 
GRB 970508 and GRB 971214 exceeds 10 52 (e.g., Waxman 1997) and 10 53 (Kulkarni et al. 1998) 
ergs, respectively, if the emission is unbeamed. The time profiles of the X-ray afterglow light 
curves are generally well fit by power laws (e.g., Piro et al. 1998; Feroci et al. 1998). 

These observations provide support for the fireball/blast-wave model of GRBs (Rees & 
Meszaros 1992). In this model, a large quantity of energy released in a small volume produces a 
GRB when the baryon loading of the fireball yields a blast wave with Lorentz factor T ~ 100-300. 
Such values of T are required to avoid significant attenuation of gamma-rays through pair 
production processes (e.g., Baring & Harding 1997). Fireballs with larger and smaller baryon 
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loading are difficult to detect because of design limitations of past telescopes (Dernier, Chiang, & 
BSttcher 1998). 

The success of the relativistic blast-wave model is due largely to the relative simplicity with 
which it explains the temporal dependence and intensity of the long wavelength GRB afterglows 
(e.g., Wijers, Rees, & Meszaros 1997; Waxman 1997; Vietri 1997; Chiang & Dermer 1998). 
Power-law decays in the afterglow light curves result from the deceleration of the blast wave as it 
becomes energized by sweeping up material from the circumburst medium (CBM). This is termed 
the external shock model of GRBs, and was originally introduced to provide a mechanism for 
converting the energy of the blast wave into radiation during the prompt gamma-ray luminous 
phase of GRBs (Meszaros & Rees 1993). 

Recent research (Fenimore, Madras, k, Nayakshin 1996, Sari & Piran 1997; Fenimore et 
al. 1998) indicates that the external shock model faces difficulties in explaining short timescale 
variability (STV) in GRB light curves and afterglows as a consequence of special relativistic effects 
resulting from the curvature of the blast wave shell. In this Letter, we construct a computer model 
to examine requirements of an external shock model to produce STV. In Section 2, we outline the 
temporal spreading problem due to the curvature of a blast wave. Our numerical model used to 
simulate this system is described in Section 3, and results are presented in Section 4. If STV in 
GRB light curves is due to an external shock model, then the required characteristics of such a 
model are outlined in Section 5. 



2. Variability from Localized Blast- Wave Emission Regions 

The STV problem in the external shock model is illustrated in Figure 1. A blast wave, 
expanding with Lorentz factor T, sweeps up material from a density inhomogeneity, or "cloud," 
located at a distance R from the explosion site. We assume for simplicity that the cloud is 
spherical with radius r. Let 9 represent the angle between the directions from the GRB explosion 
site to an observer and to the center of the cloud. The blast wave first interacts with the cloud 
at point 1, and stops interacting when it reaches point 4. The largest angular extent of the blast 
wave energized by the cloud is defined by the locations of points 2 and 3. The observer will record 
photon arrival times from each of the four points. The time delay between photons emitted from 
points 1 and 4 is 

At r = -|(i-/w)(i + z) = 2r ^ v Z) =W^(i + fl 2 r 2 )], (l) 

where the Doppler factor T> = [T(l — /3cos#)] -1 , (3 = (1 — 1/T 2 ) 1 / 2 , and z is the redshift. The 
quantity tdur = -Rmax(l + z )/(T 2 c) represents the characteristic duration of emission observed from 
a blast wave which radiates as it passes through a medium of maximum extent i? max > R- The 
expression on the right side of eq. (1) holds in the limit r>l, implying that 6^1 because most 
of the radiation is observed from a cone of angular extent^ 1/T due to the Doppler effect. For 
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constant T, the time delay between photons emitted from points 2 and 3 is 

2E r 
At a = — sm9 sin X (1 + z) = t dnr [-(2T 2 9)} , (2) 
c H 

where we assume that r/iJ<l and T » 1 in the right-most expression of eq. (2). 

The term At r is a radial spreading time scale related to the observed interval over which the 
blast wave passes through the radial extent of the cloud, and the term At a is an angular spreading 
time scale related to the angular extent of the cloud (see Fenimore et al. 1996; Sari & Piran 1997). 
The shortest temporal duration from the blast-wave/cloud interaction is defined by the maximum 
of these two time scales. Note, however, that the pulse width from an interaction event could be 
longer then the kinematic minimum time scale derived here if the blast wave does not radiate its 
swept-up energy rapidly, or if the blast- wave shell has thickness £ ;> r. In the limit that the blast 
wave is approximated by expanding collisionless particles, £ ~ R/(2T 2 ) (Meszaros, Laguna, & 
Rees 1993). The shell would be much thinner, however, if treated as a radiative hydrodynamical 
fluid that is compressed by the pressure of the external medium (Blandford & McKee 1976), and 
we assume that this is the case here. 

As can be seen by comparing eqs.(l) and (2), At a is much longer than At r except within a 
very narrow cone of extent 9 £ l/(2r 2 ). If the blast wave is uncollimated, much of the observed 
radiation is produced from interactions near 9 ~ 1/T where At a /At r ~ T. To produce STV 
with At a < 5t <C tdur therefore requires very small clouds with r/(R/T) <^ (St/t^ ur ). GRB 
time profiles commonly show pulses with dt/t^nr ^ 0.01 — 0.1; thus very small clouds which 
individually have a small surface covering factor are required for STV. The average covering 
factor is the ratio of the total area of the clouds to the area of the beaming cone, and is given 
by ~ / c 7rr 2 /[7r(i? max /r) 2 ] = I c f 2 , where I c is the number of clouds in the beaming cone and 
/ = r/(i? max /T). The efficiency for extracting energy from a blast wave is therefore roughly equal 
to 

C = min(iV c /iV C(i , 1) I c f 2 . (3) 

Here N c is the column density of a cloud, and N c d is the column density that a cloud must have 
to extract a significant fraction of the available energy in the portion of the blast wave which 
interacts with the cloud. Here we implicitly assume that a large fraction of the blast wave energy 
is converted to nonthermal electrons and then to radiation. If this is not the case, then the 
efficiency is reduced and the energy requirements increased accordingly. If the explosion energy 
Eq = 10 54 £ , 54 ergs is emitted uniformly in all directions from the explosion site, then strong 
deceleration occurs when the swept-up relativistic cloud mass Tm c equals the baryonic mass in 
the region that interacts with the cloud, where m c = 4irN c m p r 2 /3 is the cloud mass. This implies 
N cd ^ 3E /(16irr 2 m p c 2 R 2 ) = 4.4 x W w E 54 /Tl 00 R 2 7 cm" 2 , where T = 300r 300 and R = W 17 R 17 
cm. 

In order to extract the blast wave's energy efficiently by interacting with clouds, we see from 
eq. (3) that a large number of clouds with thick columns (i.e, iV c ^ N c d) are required to offset 
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their small size. If, for example, / = 0.01, then ^ 1000 clouds with thick columns within the 
Doppler beaming cone will give an efficiency ( ;> 0.1. It has been argued that pulses from the 
individual clouds would overlap (Sari &: Piran 1997; Sumner & Fenimore 1998) and fail to produce 
light curves with large amplitude variability when the efficiencies are large. We have constructed 
a model to assess this claim quantitatively. 



3. Description of the Model 

We use Monte Carlo methods to simulate the interaction of a spherically expanding blast wave 
with clouds in the CBM. The clouds are randomly distributed in angle, azimuth, and location. The 
radial distribution can be generally chosen but, for purposes of illustration, we consider a uniform 
distribution in the range R < R max = 10 17 cm. The mass and radius distributions of clouds can 
also be generally chosen, but here we consider for clarity the case where all clouds have equal 
radii. In the results shown here, we assume that the clouds have thick columns. This means that 
N c ^> 4 x 10 18 cm -2 for clouds with R > 10 16 cm; closer clouds require larger column densities, but 
the number of such clouds is extremely small. We employ a blocking factor to account for occulted 
clouds. After the locations and radii of the clouds are selected, the program determines which 
clouds have been completely occulted and eliminates them. When partial occultation occurs, we 
weight the contribution of that cloud based on the fraction of the area that is occulted. 

We use the following expression to represent the angle-dependent flux density 5(ergs cm" 2 
s _1 e~ b 1 s ) measured at observer time t ^ s from that portion of the curved blast wave surface which 
is energized by its interaction with a cloud: 

S(e ohs , Us; lUs) = K P3(1 /; )J(£) 4(^) VV^WV-s , (4) 

where 

f 1 and At = At r , if 9 < 1/{2T 2 )- 

K=\ (5) 

{ r{RTV sine sin X )~ l = (T9V)' 1 and At = At a , if (9 > l/(2r 2 ) 
(see Dermer, Sturner, & Schlickeiser 1997). In eq. (4), We approximate the pulse duration by a 
Gaussian function with a = At/(4^/2 In 2), which gives a FWHM pulse duration =At/2. In eq. 
(5), e Q bs = huohs/rrieC 2 = Ve/(l + z) is the observed dimensionless photon energy and e is the 
photon energy measured in the comoving frame of the blast wave, t c = (1 + z)R(l — [3 cos 9) / \(5c) 
when At r > At a , and t c = (1 + z)R{\ — [3 cos 9 cos x) I (P c ) when Ai r < At a The quantity J(e) 
represents the spectral power in the comoving frame. The above result is obtained by noting 
that in the limit r <^ R/T, the beaming factor for the time-integrated fluence from the radiating 
portion of the blast wave energized by the interaction is the same as that for a spherical ball of 
plasma which radiates for a comoving period of time equal to ti — t\ = 2r/j3Tc. 

We approximate the comoving radiant emission by a power law such that J(e) = Joe _Q in 
the energy range ei < e < €2, where a is the photon energy index. To set the normalization 
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for Jo, we note that due to the increasing surface area of the expanding blast wave, the 
maximum energy that could be radiated as a result of the blast wave/cloud interaction is just 
E c = [Eqt 2 / (AR 2 )] mm(N c /N cc i, 1). We assume that this energy is radiated promptly, that is, the 
interaction takes place in the radiative regime. The relationship between the measured energy 
fluence Fe (ergs cm~ 2 ) and the total energy emitted is Fe = (1 + z)E c /(4ird 2 L ), where Jl is the 
luminosity distance. From this relation, we obtain the result Jo = f a [E c /T(t2 — = f a (cE c /2r), 
where / Q = (1 - a)/(e£- a - ej- a ). 



4. Results 

In our calculations, we choose parameter values motivated by the observations of the candidate 
host galaxy of GRB 971214 at redshift z = 3.42 (Kulkarni et al. 1998). Its total energy release, 
including prompt radiation and afterglow, could exceed 10 54 ergs if the burst source is unbeamed. 
Thus we choose a value of dE/d£l = 10 53 ergs sr _1 for our standard GRB energy release, which 
implies a value of E54 = 1.3 if the burst energy is isotropically radiated. For the calculations 
shown here, we let T3oo = 1, z = I (which implies Jl = 1.67 x 10 28 cm for a Hubble constant of 
65 km s _1 Mpc -1 and qo = 1/2), a = 1, €\ = 1CP 4 , and 62 = W 2 . In order to conserve computing 
time, all clouds are placed at angles 9 < kg/T with kg = 3, as we find that clouds located at larger 
angles make a negligible contribution to the observed flux due to the strong Doppler beaming. 
The clouds are located in a shell between -R m i n = 10 16 and i? max = 10 17 cm. 

Fig. 2 shows results of our calculations for a range of different cloud sizes cit £obs — 

1. The 

number of clouds within the Doppler beaming cone 6 < 1/T is chosen in order to obtain a 10% 
efficiency as defined by eq. (3). In curves (a), (b), and (c) we let / = r/(R max /T) = 0.1,0.03, 
and 0.01, implying a total number of clouds kgl c = 9I C = 90, 1000, and 9000 in the 
respective calculations. From the expression for N c d following eq. (3), this implies cloud 
densities n c ;> 100£ , 54/(r3oo-R'i7/) cm~ 3 . By examining the curves, one can verify that the 
calculated energy fluences agree with the expected energy fluence = 0.1(1 + z)4tt x 10 53 ergs 
si~ l /[4ird\ ln(e2/ei)] ~ 5 x 10~ 6 ergs cm~ 2 for a 10% efficiency. 

From eq. (2), we see that a typical pulse duration is ~ At a ~ 2r(l + z)/(Tc) when 9 ~ 1/T. 
The cloud size in curve (a) has a value of r = 3.3 x 10 13 cm, implying St = 15 s when 9 ~ 1/T, 
comparable to the value of tdur = 74 s. Thus only limited structure in the light curves are 
seen. For the cases with smaller clouds shown in curves (b) and (c), the light curves become 
increasingly variable. In many but not all cases, a very narrow pulse with large amplitude is 
obtained as a result of a cloud which happens to be located within an angle 9 <J 1/r 2 of the 
observer's line-of-sight. These early bright peaks might not be so apparent if the blast-wave 
thickness also determines the pulse duration. The light curves are most variable during the period 
~ (1 + z)(R mSLX — R m - m ) I (2r 2 c) ~ 33 s. The tendency for longer pulses to be found in the simulated 
light curves at later times is an artifact of the bounded cloud ranges in the numerical simulation. 
This effect is ameliorated when more complicated cloud profiles without sharp boundaries are 
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considered. All these simulated GRBs would be detectable with BATSE on CGRO out to redshifts 
z 4, noting that BATSE triggers on a vF v flux of ^ 10 -7 ergs cm~ 2 s -1 . Higher efficiency events 
could be detected from sources at even larger redshifts (see inset to Fig. 2). 

The mean number and FWHM durations of pulses for an ensemble of GRBs corresponding 
to cases (a), (b) and (c) in Fig. 2 are 7.9, 19 and 42, and 4.8 s, 1.5 s and 0.42 s, respectively. We 
consider only those pulses which either rose from or decayed beyond 50% of their peak value, and 
used twice the half- maximum duration when only one of the aforementioned cases was satisfied. 
Norris et al. (1996) has decomposed 41 bright GRB light curves to find a distribution of pulse 
durations in the range 0.2-2 s. As shown in the right inset to Fig. 2, cloud sizes intermediate to 
cases (b) and (c) are in accord with the distribution deduced by Norris et al. A distribution of 
cloud sizes and the addition of background noise in the simulations could improve the comparison. 

5. Discussion and Summary 

In the external shock model considered here, STV results from the energy of a blast wave 
being efficiently extracted by sweeping up matter from the CBM in the form of small clouds 
(compare the approach of Shaviv h Dar 1995). If the clouds are distributed uniformly around the 
burst source, then a total cloud mass of 

M c 4r 2 / c • ^m p N cMn r 2 £ 3 x lO" 4 ^ ( f 54 ) M (6) 

6 1 300 ^16 

is required. Here the minimum cloud column density Ar c min is derived for clouds at 10 16 i?i6 cm. 
In contrast, the minimum mass of the surrounding intercloud medium (ICM) required to extract 
a significant fraction of the blast wave's energy is Eq/{T c 2 ) «6x 10 _6 £ , 54/r 2 00 M & . More mass 
must be in the form of small clouds to ensure that their columns are thick at all values of R, but 
the total mass is remarkably small in either case. The contrast between the ICM density uicm 
and the cloud density n c required to obtain STV is njc]vi/n c <C (4/£/T) <C 1. Whether such 
clouds can be stably confined or must be continuously reformed is unclear. 

We consider arguments presented against the external shock model (Fenimore, Ramirez, & 
Sumner 1998; Fenimore et al. 1996; Sari & Piran 1997; Sumner & Fenimore 1998). The strongest 
objection is that the overlapping light curves will "wash out" large amplitude variabilty. According 
to eq. (2), most clouds produce pulses of duration 5t ;> 2ft^ LVir {T0). Because most clouds in the 
Doppler cone are located at 9 ~ 1/T, STV only occurs when St/t^ UT ~/<Cl. To obtain a given 
efficiency, eq. (3) shows that I c = Q/ f 2 clouds with thick columns are required. The number of 
overlapping clouds in the time element St is « IciSt/t^) ~ fl c ~ (/f. This would apparently 
yield statistical fluctuations of order (f /C) 1 ^ 2 which, for / = 0.01 and ( = 0.1, is at the 30% level 
rather than at the factor-of-2 level often seen in GRB light curves. 

This argument overlooks the fact that the small fraction of clouds at 6 <C 1/r makes a 
disproportionate contribution to the amplitude and variability of GRB light curves. Only 1% of 
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all clouds in the Doppler beaming cone lie in the range < 9 < l/(10r), yet the duration of pulses 
from these clouds is shorter by an order-of-magnitude than the duration of pulses from clouds at 
9 = 1/r. Moreover, the amplitude of the pulses varies according to the factor ~ p( 2 + Q ) /(T9), 
which represents a factor « 80 for clouds at 6 < l/(10r) when compared to clouds at 9 = 1/r. We 
therefore see that 1% of the clouds produce « 8% of the total fluence in very short duration, very 
large amplitude pulses. The amplitude enhancements continue to increase oc 1/9 until 9 <^ 1/(2T 2 ), 
or until the pulse width becomes limited by non-kinematic effects. If clouds are found in layers, 
then the enhanced contribution of on-axis clouds might also overcome the difficulty in producing 
gaps in GRB light curves (Fenimore, Ramirez, & Sumner 1998). 

The efficiency £ defined in eq. (3) represents the most optimistic case where all the blast wave 
energy in the region energized by the cloud is transformed into radiation. Because the blast wave 
decelerates strongly in the region of interaction, STV can be produced through Doppler deboosting 
even if the energized electrons do not lose most of their energy through synchrotron processes 
(Chiang 1999a, b). In this case, of course, the efficiency is much less than obtained through eq. (3). 

In summary, we have examined a model in accord with the third scenario described by 
Fenimore et al. (1996), where STV in GRB light curves is produced by inhomogeneties in the 
CBM. This breaks the condition of local spherical symmetry, shown by Fenimore et al. to produce 
light curves inconsistent with observations, but requires the presence of many small clouds with 
thick columns. As shown here, STV in GRB light curves can be produced by the external shock 
model under such conditions. Further comparisons of GRB data with model light curves will be 
necessary, however, to establish whether GRB light curves are tomographic images of the density 
distributions of the medium surrounding the sources of GRBs 
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Fig. 1.- 



Geometry of the interaction between a blast wave and a cloud in the circumburst medium. 
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Fig. 2. — Simulated light curves from blast- wave/cloud interactions. The cloud radii range from 
3.3 x 10 13 cm in curve (a) to 3.3 x 10 12 cm in curve (c), and the total number of clouds is estimated 
from eq. (3) in order to extract 10% of the blast wave energy (see text). Successive curves are offset 
by 50 s each for clarity. Left inset shows a calculation for 3.3 x 10 12 cm clouds giving ~ 100% 
efficiency. Right inset shows the distributions of pulse widths, normalized to unity, for curves (case 
b; dashed histogram) and (case c; solid histogram). The early narrow peak in graph (c) extends to 
3 x 10~ 6 ergs cm -2 s _1 , but is covered by the inset. 



